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Abstract. We present a minimal continuum model of strongly adhering cells as active 
contractile isotropic media and use the model to study the effect of the geometry 
of the adhesion patch in controlling the spatial distribution of traction and cellular 
stresses. Activity is introduced as a contractile, hence negative, spatially homogeneous 
contribution to the pressure. The model shows that patterning of adhesion regions can 
be used to control traction stress distribution and yields several results consistent with 
experimental observations. Specifically, the cell spread area is found to increase with 
substrate stiffness and an analytic expression for the dependence is obtained for circular 
cells. The correlation between the magnitude of traction stresses and cell boundary 
curvature is also demonstrated and analyzed. 
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1. Introduction 

Living cells actively sense and respond to the physical geometry and stiffness of their 
environment, which in turn affects a variety of cellular processes, such as growth, 
differentiation, morphogenesis, spreading and motility [Ij. Cell-matrix adhesion is 
mediated by integrin complexes, referred to as focal adhesions, that bind to specific 
ligands on the underlying matrix. Focal adhesions are mechanically linked to the 
actomyosin cytoskeleton inside the cell, that in turn generates contractile forces on the 
extracellular matrix. The interplay between substrate stiffness, intracellular contractility 
and the extracellular adhesion forces controls the cell morphology and its mechanical 
behavior. For instance, cells adhering to soft substrate are generally found to spread 
less and have round morphology, while cells on stiff substrates have greater spread area 
with more branched shapes [2]. 

Powerful techniques have been developed in recent years to measure the traction 
forces exerted by adherent cells on compliant substrates [3] . Traction Force Microscopy 
is used to probe the traction stresses exerted by cells on continuous elastic gels. 
The stresses are inferred from measurements of the displacements of fiducial markers 
embedded in the gel before and after cell detachment [4l|5]. In a second technique cells 
plated on microfabricated pillar arrays induce bending of the elastic micropillars. The 
traction forces are then obtained from by assuming a linear Hooke's law relation between 
the measured bending and the forces [6] . These experiments have demonstrated that the 
mechanical response of adherent cells is controlled by a complex interplay of substrate 
stiffness and geometry, myosin activity and extracellular matrix proteins. Adhesive 
micro patterning has also been used as a tool for both controlling cell shape and study 
the interplay between shape and cytoskeletal organization and architecture [7J. These 
studies have shown that when strongly adhesive patterns force the cell boundary to 
exhibit regions of high curvature, traction stresses tend to be concentrated in these 
regions, while stress fibers develop along cell boundaries linking non-adhesive zones, 
confirming the crucial role of the cytoskeletal contractility and architecture in controlling 
cellular stresses and morphology [8]. 

Network models of the contractile cytoskeleton have been used to describe the 
relation between force distribution and shape of adherent cells and are discussed in a 
recent work by Guthardt Torres and coworkers [9j. These models include networks of 
Hookean springs as well as cable networks that incorporate the asymmetry of the elastic 
response of biopolymers such as filamentary actin to compression version extension, 
with and without the explicit inclusion of contractility. In particular, the active 
cable network introduced in [9j reproduces the arc morphology of cell boundaries 
pinned by strong local adhesions that has been seen in experiments. In this paper 
we consider, in contrast a continuum model of cells as linear, active elastic media 
and demonstrate that the introduction of activity simply as a spatially homogeneous 
contractile, hence negative, contribution to the pressure is sufficient to reproduce 
the spatial inhomogeneous distribution of traction and cellular stresses observed in 
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experiments for a number of cell geometries. An interesting extension of our work 
will be to introduce nonlinearity in the continuum model to describe an asymmetric 
response to compression and stretching. This asymmetry, arising from the nonlinear 
force-extension curve of actin filaments, is known to be important in controlling the 
contractile behavior of isotropic gels [TOl [H] and may alter the stress distribution in 
adhering cells. 

In the next section we introduce our continuum model of adherent cells as active 
contractile elastic media. We then use the model to study the effect of the geometry 
of the adhesion region on controlling the spatial distribution of stresses in the cell. 
The model can be solved analytically for a circular cell, where we obtain an expression 
for the cell spread area as a function of substrate stiffness and show that our results 
compares favorably to experiments (inset of Fig. [2]). The cases of elliptical, square and 
triangular cells are solved numerically. We show that the geometry of the adhesive 
region strongly affect the stress distribution, with traction stresses concentrated in 
regions of high curvatures or at sharp corners (Fig. [T]). In section 3.3 we provide an 
analytical argument that quantifies the correlation between traction stress magnitude 
and curvature of the cell boundary and discuss in section 3.4 the relative roles of shear 
and compressional deformations in controlling the stress distribution. We conclude with 
a brief discussion. 



2. Adherent cell as a contractile gel 

We consider a stationary cell adhering to an elastic substrate via stable focal adhesion 
complexes. We further assume that the cell has attained its optimum spread area on 
the substrate, with an average height h much thinner than its perimeter. In mechanical 
equilibrium, the condition of local force-balance translates to df^aa/s = 0, where a is 
the three-dimensional stress tensor of the cell with greek indices taking values x, y and 
z. For a thin cellular film we average the cellular force-balance equation over the cell 
thickness h. In-plane force balance is given by 

djaij + d^aiz = , (1) 

with i^j denoting in-plane coordinates. We assume that the top surface of the cell is 
stress free, criz{rj_^z = /i) = 0, whereas at the cell-substrate interface z = 0, the cell 
experiences lateral traction stresses given by cr^z(r^,z = 0) = Yui{r_i^z = 0). Here, 
Y denotes the substrate rigidity parameter, representing the cell-substrate anchoring 
strength, and u(r^,z) is the in-plane deformation field of the cellular medium. The 
thickness- averaged force balance equation then reads fi2\ [T3] . 

hdjaij = Yui , (2) 

where a^j(r^) = J^{dz/h)aij{r_\_^ z). In the following we will drop the overbear 
indicating the average and refer to thickness averaged quantities throughout. The 
quantity = Yui is a stress in three dimensions, i.e., a force per unit area. It 
describes the in-plane traction force per unit area that the cell exerts on the substrate. 
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The assumption of local elastic interactions with the substrate strictly holds on elastic 
substrates that are much thinner than the lateral size of the cell [I4j or on micropillar 
substrates [13j. The substrate rigidity parameter Y depends on the stiffness of the 
underlying substrate as well as on the density pf and stiffness fc/ of focal adhesions. For 
an elastic substrate of shear modulus /x^ and thickness Y takes the simple form [T4], 

kfpf Ps/hs ' 

We model the cell as an isotropic and homogeneous elastic material with additional 
internal active stresses due to actomyosin contractility. The constitutive relation for the 
cellular stress tensor is then given by, 

V • u 5ij + djUi + diU^ + Ga^ij , (3) 

where E and v denote the Young modulus and Poisson ratio of the cellular material, 
respectively. Actomyosin contractility is modeled as a negative contribution to the 
local pressure, corresponding to > 0. The assumption of linear elasticity is valid 
on time scales shorter than cytoskeletal turnovers, that are indeed slowed down by 
strong adhesion to the substrate. Equations ^ and ([3]), subject to the boundary 
condition cr^j-rijl^^ = 0, wtih Vt the cell boundary and n the outward unit normal on 
ri, completely describe the equilibrium of an adherent cell. As a consequence of the 
stress free condition at the lateral cell boundary, the net traction force transmitted by 
the cell to the substrate vanishes, i.e., /^rf^r Yui = §^ds crijUj = 0. It is instructive 
to consider two limiting cases for the anchoring strength. When the cell is rigidly 
anchored onto the substrate, corresponding to y ^ oc, we find u = 0, defining the 
reference state for elastic deformations. In our model the reference cell shape is then 
dictated by the geometry of the adhesion patch, which can be controlled in experiments 
by micropatterning substrates by adhesion proteins. In contrast, when y ^ 0, the cell 
does not adhere to the substrate and the equilibrium state is uniformly contracted state, 
with a density enhancement 5p = — V-u = aa{l + iy){l — 2iy)/ E{l — iy). In the following, 
we investigate analytically and numerically solutions of the cell elasticity equations ^ 
and ([3]) subject to stress-free boundary conditions in various planar geometries. 

3. Results 

3.1. Spatial distribution of traction stresses is sensitive to adhesion geometry 

The spatial distribution of traction stresses exerted by cells on substrate and the 
corresponding organization of stress and deformation inside the cell are affected by 
the geometry of adhesive patterns. Using micropatterning techniques, cell shapes can 
be constrained to adhere to controlled geometrical patterns [151 US] • In our model the 
shape determined by the pattern in the limit of infinite adhesion strength provides the 
reference shape for the cell. Here we investigate four reference cell shapes: circle, ellipse, 
square and equilateral triangle. These are chosen to have the same reference area but 
different perimeters. The case of a circular cell can be treated analytically, as described 
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Figure 1. Equilibrium cell shapes for various adhesion patterns : Circle (top left), 
ellipse (top right), square (bottom left) and equilateral triangle (bottom right). The 
color map indicates magnitude of the traction |T| = y|u|, and the arrows demote the 
direction of the traction vectors. The reference shapes for all the four patterns have an 
equal area of 1000 jam^. The other parameters are: E = 1 kPa, u = 0.4, da = I kPa, 
fis = 10 kPa, hs = 30 //mm, h = 0.2 jam. 

below. For the other shapes the elasticity equations ^ and ^ are solved numerically 
using the MATLAB pde toolbox. We assume the contractility to be uniform and of 
order of the cellular Young's modulus. Heatmap of traction stresses are shown in Fig. [TJ 
In all cases the traction stresses are concentrated at the cell periphery, irrespective of 
the reference shape. The magnitude of the local traction stress is, however, higher in 
regions of high curvatures or at sharp corners. 

For a circular cell, Eqs. ^ and ^ can be solved analytically [I3lll7]. Assuming 
in-plane rotational symmetry, it is convenient to use polar coordinates r and 0, denoting 
radial and angular coordinates, and demand that no quantity depend on 9. The equation 
for the radial displacement Ur about a circular reference state of radius i?o, is then given 

by 

r^d'^Ur + rdrUr - (1 + r^/il)ur = , (4) 
where the penetration length ip describes the localization of traction stresses at the cell 
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boundary. It is given by : 

Eh{l 



f = 



and is essentially controlled by the ratio of the cell stiffness ^ to the substrate rigidity 
^ Y . The penetration length is short on stiff substrates and increases with decreasing 
substrate rigidity. The solution of Eq. Q with the boundary conditions a^^(r = i?o) = 
and 2i^(r = 0) = is given in terms of modified Bessel functions of the first kind as, 

'{l + i^){l-2u) 



Ur{r) 

with g{s) = [sIo{s) 



l-2u 



E{1 



h{r/e,)g{Ro/i 



pj 



(6) 



As anticipated, the deformation Ur vanishes for all 



r when y ^ oc, when the adhering circular cell is maximally spread and has its largest 
undeformed radius Rq. 



3.2. Cell spread area is sensitive to substrate stiffness and active contractility 

The optimal spread area of the cell is controlled by the interplay between cell 
contractility, as described by the active pressure a^, and the traction forces on the 
substrate. In the case of a circular cell, where the deformation induced by adhesion is 
given by Eq. ([g]), the steady state cell area is given by, 

A = 7r{Ro + u{Ro))\ (7) 

with Rq the reference radius corresponding to the maximal spread area — 
attained on an infinitely rigid substrate, where Ur{r) = 0. To make contact with 
experiments, we investigate the ratio A/A^o^ the relative cell spread area, as a 
function of substrate stiffness and contractility. On stiff substrates, where Rq ^ £p^ 
i.e., the traction stress extends over a length much smaller than the reference cell 
radius, Ur{Ro) — —cfa^p/B^ where the compressional modulus B is given hy B = 
E{1 — u)/ [{1 + — 2iy)]. The relative spread area then takes the simple form 

A,A^ ^ (l _ Lettmg Y ^ ,./h., wa note that increa.,„g substrate 

stiffness increases relative spread area, with A/A^^ ^ 1 as /x^ ^ oc, in qualitative 
agreement with experiments [21 [181 [19]. contrast, increasing the contractile pressure 
aa reduces the optimal cell spread area, consistent with the experimental observation 
that myosin-II activity retards cell spreading [20j. To make a quantitative comparison 
with experiments, we fit Eq. ([7| to experimentally reported data on the projected area of 
cardiac myocytes cultured on N-cadherin coated PA gels of varying stiffnesses [l9] . Here 
the maximal spread area is taken to be equal to the cell projected area on a glass 
substrate (shear modulus ^ 30 GPa), which is 690/xm^. The fit, shown in the inset of 
Fig. [2} is obtained using the active contractility and the cellular Young's modulus E 
as the fitting parameters. A least-square fit gives us = 911 Pa and a a = 1589 Pa. This 
value is consistent with the value used in Ref. [ll7j to fit the measured value of the surface 
tension of a colony of epithelial cells. Figure [2] also demonstrates the competing roles 
of contractility and adhesion in controlling optimal cell shapes for a chosen triangular 
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Figure 2. Optimal shape of a triangular cell for different values of the active pressure 
aa and the substrate shear modulus /ig, with E = 1 Pa. The color map represents the 
magnitude of the displacement vector |u| (proportional to the traction force) about an 
equilateral triangular reference shape of area 1000 /im^. The cell spread area increases 
with increasing substrate stiffness and decreases with increasing aa- Inset : Least- 
square fit of the relative cell spread area A/Aoo obtained from our theory (see Eq. ^) 
to the experimental data reported in Ref. [19]. The fitting parameters are E = 911 Pa 
and aa = 1589 Pa. The other parameters are: v = 0.4, hs = 30 /im, h = 0.2 /im. 



reference state. On softer substrates the triangular cell retains its topology and contracts 
by an amount proportional to a^, whereas on stiffer substrates the corners tend to form 
protrusions. 

3.3. Traction stress magnitude correlates with cell boundary curvature 

When the boundaries of the adhesion pattern exhibits non-uniform curvature, the 
traction stresses are higher at regions of high curvatures. This is seen for example 
in Fig. [T] for the case of an elliptical reference shape. To justify this we propose a simple 
analytical argument based on local force balance. Consider a thin slice of cellular 
material at the cell periphery of width comparable to penetration length £p and arc 
length RAO much less than the cell perimeter (Fig. |3](a)), with 1/R the local curvature 
of the cell element. At the outer edge of this element, the only force on the cell is the 
reaction to the traction by the cell on the substrate traction, of areal density — T, with 
T = Yu. This yields an outward total force on the outer edge of the cell element of 
magnitude TRAOip^ with T > 0. At the interior edge, the cellular element experiences 
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(a) T(traction) 




Figure 3. (a) Force-balance on a thin slice of cellular material at the cell boundary, 
(b) Force-balance at a generic sharp corner with opening angle (j). (c) Traction stress 
magnitude at the cell edge as a function of the local curvature n for the elliptical cell 
of Fig.[l] 



a contractile force of magnitude cr^(i? — £p)A9£p^ where is the normal stress pulling 
the inner contour inwards and has contributions from active as well as passive elastic 
stresses. The lateral stresses at contributes to an effective line tension atipRAO of the 
cell element. Due to the curvature of the boundary element, the line tension generates 
an inward Laplace pressure of magnitude at£p/R. Local balance of forces then yields, 

TRAOip -(Jn{R- QAOip = RAeipGt^ . (8) 

R 

The above law can be written down in a compact form as, 

T = Gn + [Cft- Crn)£pl^ , (9) 

with tv = 1/R^ the local curvature of the boundary element. Equation ^ then tells us 
that local magnitude of traction increases linearly with increasing boundary curvature. 
The lateral and normal stresses at and can be expressed in terms of the local cellular 
stresses in polar coordinates as at = aee — dearo and a^ = cfrr- The linear dependence 
of T on strictly holds in the limit £p/<i <C 1. In addition, non-local elastic interactions 
can also affect the dependence of traction magnitude on local curvature. Figure. [3](c) 
shows the dependence of the magnitude of the traction stress at the cell boundary on 
local curvature for an elliptical cell as shown in Fig. [T} For low the traction stress 
magnitude increases linearly with k before reaching a plateau at higher values of 
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Figure 4. Cell shape anisotropy correlates with internal stress anisotropy. (a) 
Heatmap of internal compressive stress a (left) and maximum shear stress ag (right) 
corresponding to various reference shapes : circle, ellipse, square and equilateral 
triangle. The reference shapes all have an equal area of 1000 /im^. (b) Average 
maximum shear as a function of eccentricity e for elliptical cells of same reference 
area (1000 jim?). Equilibrium shapes with colorplot of ji are given as plot markers. 
Parameters : E = 1 kPa, v = 0.4, (7^ = 1 kPa, /i^ = 10 kPa, hg = 30 /im, h = 0.2 jam. 



When the cell boundary exhibits a sharp corner with opening angle 0, as shown in 
Fig. [3](b), the local force-balance is given by, 

T = a, + 2a, cos (0/2) , (10) 

where acts along the bisecting line of the corner. Hence smaller the opening angle, 
the larger is the traction force. 

3.4' Internal stress anisotropics induced by anisotropic adhesion geometry 

The spatial distribution of internal stresses aij within the cell depends on cell shape, 
which is in turn controlled by the geometry of the adhesive region. Experimentally 
a{x^y) can be obtained from the measured distribution of traction stresses T(x,^), 
inverting the local force-balance condition djaij = Ti [21]. The elasticity equations 
Eqs. ^ and ^ can be recast as a single partial differential equation for the internal 



Controlling cell-matrix traction forces by extracellular geometry 



10 



stress tensor a^j, given by 

[didkCTkjf + Sijaa = cFij + —-^y—^ij i^kk " 2a«) , (H) 

where [...]^ denotes symmetrization with respect to indices that are not summed over, 
i.e., [didko'kj]^ = I [didk^kj + djdk^ki]- We have investigated numericaUy the solution of 
Eq. (11) with stress free boundary condition crijUj = 0. To understand the role of shear 
and compressional deformations in different geometries, it is instructive to diagonalize 
the stress tensor and display the results in terms of linear combinations of the eigenvalues 
Gi and (72. The sum a = \{cri + 72) is simply the trace of the stress tensor and describes 

compressional deformations. The difference = — 72! = ^Jl^xx ^yyY + ^^ly-) 
controlled by normal stress a^x — (^yy and shear stress For an isotropic reference 
shape, such as the circle, ai = 02 and Gs = 0, whereas for anisotropic shapes such as 
the ellipse, one expects nonzero values for the local maximum shear a^. 

Fig.|4](a) shows heatmaps of the spatial distribution of a and for various reference 
shapes - circle, ellipse, square and equilateral triangle. Irrespective of the shape of the 
adhesion geometry, a is maximum at the cell center, indicating build-up of compressive 
stresses. The compressional stress a always vanishes at the boundary, and it does so 
more rapidly at regions of high curvature or at sharp corners. In contrast, the shear 
stress Gg is identically zero for isotropic shapes, defined as those that have a gyration 
tensor that is diagonal, with equal eigenvalues. The circle, triangle and square are all 
in this class. Local stress anisotropy as measured by Gg is nonzero for elliptical shapes 
and shear stresses build up at the center of the ellipse. The shape anisotropy of ellipses 
can be quantified by their eccentricity e = ^1 — with a and h the semi-major 

and semi-minor axes. Figure |4](b) shows the spatial average of Gg over the area A of the 
cell, defined as = ^ Ia^^^ cr^, as a function of the eccentricity e. The average shear 
stress Gg increases with e with a sharp rise as e ^ 1, indicating a positive correlation 
between geometrical and mechanical anisotropy in adherent cells. 



4. Discussion 



We have used a continuum model of an adherent cell on a substrate as an active 
contractile medium to study the role of adhesion geometry in controlling cell shape, 
cell spreading and the spatial distribution of traction stresses. The central input of the 
model is the cell contractility or activity a^, a negative contribution to the pressure 
that enters the constitutive equation for the cellular material. In general, Ga will be 
determined by the concentration and activity of myosin proteins cross linking the actin 
cortex and controlling the formation of stress fiber. In our model Ga is assumed to be a 
constant parameter, to be determined by fitting experiments. We consider cells adhering 
to flat substrates that have been patterned with adhesive patches, consisting for instance 
of fibronectin coatings, of specific geometry and examine the role of the geometry of the 
adhesive patch in controlling the spatial distribution of stresses in the cellular material. 
The reference state for our cell is the limit of infinitely strong adhesion, where the cell 
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shape and lateral extent and determined entirely by the shape and size of the adhesive 
patch. For finite adhesion strength, cell elasticity and contractility yield deviations 
form this reference state. We restrict ourselves to considering continuous or densely 
spaced adhesion sites. For discrete or sparsely distributed adhesion sites, non-adherent 
segments in the cell boundary could likely exhibit morphological transitions induced by 
contractile activity and substrate stiffness, as discussed in a recent model for an adherent 
cell as a contractile film [22]. In agreement with experimental observations, we find that 
cells spread more on stiff substrates and we provide an expression for the cell area versus 
substrate stiffness for the case of a circular cell. We show that this expression does fit 
data in cardiac myocytes (see inset of Fig. 2). Significant is the fact that the value of 
obtained from the fit is of the same order as the value used in our earlier work to fit the 
measured surface tension of adherent epithelial cell colonies [17j. We also demonstrate 
that strong traction stresses correlate with region of high cell boundary curvature, again 
in agreement with experimental observations. 

Understanding the relation between cell morphology, the cell's mechanical response 
and cell fate is an important question in cellular biophysics. Our simple model 
highlights the correlation between the geometry of adhesion sites and cell morphology 
and demonstrates that traction forces by cells can be tuned by controlling the geometry 
of adhesive regions. An important open question not addressed by this simple model 
where the adhesive patch geometry solely controls the cell shape is how cell morphology 
is determined by the interplay of cell-substrate adhesion and dynamical reorganization 
of the cytoskeletal architecture in response to the adhesion stimulus. To understand this 
it will be necessary to incorporate the dynamical feedback between actin reorganization 
and adhesion kinetics. 

This work was supported by the National Science Foundation through award DMR- 
1004789. 
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